What is claimed is: 



CLAIMS 



1. A method for determining interaction sites in a protein, 
5 said method comprising the steps of: 

identifying in said protein ionizable amino acid residues 
having anomalous predicted theoretical microscopic titration curve 
behavior; and 

clustering certain of said identified ionizable amino acid 
10 residues into putative interaction sites. 

2. The method of claim 1, wherein said identifying and 
clustering steps comprise the steps of: 

(a) obtaining a three-dimensional structure of said protein; 
15 (b) calculating an electrical potential function for said 

protein structure; 

(c) calculating a titration curve for each ionizable residue 
in said protein; 

(d) evaluating the shape of said titration curve for each 
20 said residue; 

(e) identifying any said residue with a perturbed titration 
curve in comparison with other residues of the same kind; and 

(f) identifying any said perturbed titration curve residues 
that are in a cluster, wherein the existence of said cluster 

25 indicates that the residues in said cluster are at an interaction 
site in said protein. 

3. The method of claim 2, further comprising the step of 
identifying any said perturbed titration curve residues that do 

30 not fall into any said cluster. 

4. The method of claim 2, wherein, in step (a), the three- 
dimensional structure of said protein is known. 
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5. The method of claim 2, wherein, in step (a), the three- 
dimensional structure of said protein is theoretically determined 
from the amino acid sequence of said protein. 



10 



6. The method of claim 2, wherein, in step (f), said residues 
are identified as being in a cluster by virtue of being located in 
physical proximity to each other in said three-dimensional 
structure of said protein. 

7. The method of claim 5, wherein said physical proximity has a 
distance of less than 15 A. 



8. The method of claim 5, wherein said physical proximity has a 
15 distance of less than 10 A. 

9. The method of claim 5, wherein said physical proximity has a 
distance of less than 7 A. 

20 10. The method of claim 5, wherein said physical proximity has a 
distance of approximately 6 A or less. 

11. The method of claim 2, wherein said perturbed titration 
curve is an elongated titration curve where partial protonation is 

25 predicted to extend over a wide pH range. 

12. The method of claim 2, wherein said step of identifying said 
residue with a perturbed titration curve can be performed by 
visual inspection, statistical analysis or automated 

30 classification. 

13. The method of claim 2, wherein said interaction site is 
further identified as a catalytic site. 
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14. The method of claim 2, wherein said interaction site is 
further identified as a recognition site. 

5 15. The method of claim 2, wherein said interaction site is 
further identified as a cof actor binding site. 

16. The method of claim 2^ wherein said interaction site is 
further identified as a ligand binding site. 
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